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Let  y(*)  be  a stationary  mixing  process  and  J (t)  an 
approximation  to  a random  impulsive  process.  Kurtz's  {1]°  results  on 
approximation  of  a general  semigroup  by  a Markov  semigroup  are  used  to 
prove  (weak  and  a similar  type  of)  convergence  of  the  solutions  to  (1.1) 
and  (1.2f  to  jumping  diffusions.  Previous  results  are  generalized  in 
various  ways.  The  case  of  unbounded  y(^)  is  also  treated  as  is  the 
combined  jump-diffusion  case.  Also,  a limit  theorem  for  an  integral 
with  respect  to  ^approximate  white  noise"  in  terms  of  an  I to  integral 
is  given.  The  method  has  the  advantages  of  generality  and  relative 
ease  of  use. 
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1.  Introduction 


In  (1),  Kurtz  gave  some  fairly  genaral  semigroup  methods  for 
showing  convergence  of  a sequence  of  non-Markov  process  to  a Markov 
process,  either  in  the  sense  of  weak  convergence  or  in  the  sense  of 
convergence  of  finite  dimensional  distributions.  Let  y(-)  denote 
(a  Euclidean  space  valued)  right  continuous  strong  mixing  stationary 
process.  For  each  e > 0,  define  yL(t)  = y(t/t^),  and  for  suitable 
F , G , define  the  process  xL(*)  by 


F (x 


x.  = 


t'*t 


+ G(xt,yt) 


x0  x0 
y ( t ) e r1 


e 

m ’ 


m 


(1.1) 


Khazminskii  [2],  Papanicolaou  13],  Papanicolaou  and  Kohler  [4]  and 
Blankenship  and  Papanicolaou  [5]  have  all  treated  the  problem  of  weak 
convergence  of  xl(‘)  to  a diffusion.  The  problem  is,  of  course, 
closely  related  to  the  original  problem  of  Wong  and  Zakai  (6) . In 
this  paper,  Kurtz’s  results,  (together  with  a technique  exploited 
in  references  (3)  and  (5))  will  be  used  to  get  similar  types  of 
results  under  conditions  which  are  weaker. + The  method  of  proof 
has  the  great  advantage  of  being  quite  straightforward  and  easy  to 
use  for  both  the  diffusion  and  jump-diffusion  cases. 


We  also  treat  limits  of  systems  of  the  type 


_ , e E, 
.E  = F(Xt'yt) 
Xt 


T-*-  * + £ VxtIJi,t 


t.  e 


(1.2) 


where 


f\e  „ 

J.  ds  is  an  approximation  to  a pure  jump  process,  and 
J n 1 ' s 


we  obtain  a limit  which  is  a jumping  diffusion. 


References  (2)  and  (3)  allowed  F and  G to  depend  also  on  t. 
At  the  expense  of  extra  detail,  this  case  could  be  handled  by 
our  method. 


Sections  2 and  3 recapitulate  Kurtz's  method  of  proving  con- 
vergence of  finite  dimensional  distributions  and  tightness,  resp. 

The  results  are  recapitulated  partly  for  the  sake  of  self-containment , 
and  partly  to  state  the  precise  form  in  which  they  are  to  be  used. 
Section  4 states  the  assumptions  used  in  Section  5,  which  gives  the 
result  for  limits  of  (1.1)  when  y(«)  is  bounded,  a restriction 
also  used  in  the  past  references.  Theorem  3 gives  a result  which 
is  useful  in  approximating  stochastic  integrals  with  respect  to  a 
Wiener  process  by  ordinary  integrals.  Such  results  are  needed  for 
identification  and  related  problems  (see  Balakrishnan  [8],  [9]).  A 

result  on  convergence  of  finite  dimensional  distributions  for 
unbounded  y(")  is  given  in  Section  6,  and  tightness  for  unbounded 
y(*)  is  proved  in  Section  7.  Section  8 deals  with  the  relatively 
simple  case  where  the  input  is  an  approximation  to  a random  impulse 
process,  and  (1.2),  an  approximation  to  a jumping  diffusion,  is  treated 
in  Section  9.  The  result  in  Sections  6 and  7 cover  the  much  used 
case  where  y ( - ) is  a Gaussian  diffusion. 

A method  similar  to  that  of  Section  5 is  outlined  in  Section  4, 
[3],  for  the  problem  of  showing  convergence  of  finite  dimensional 
distributions  for  the  bounded  y ( • ) case.  The  results  there  are 
not  in  a particularly  usable  form,  and  actually  require  more 
smoothness  of  F and  G than  needed  here  since  partial  differential 
equation  methods  are  ultimately  used  there.  Here,  we  do  not  need  to 
solve  or  even  to  approximate  solutions  of  partial  differential 
equations,  but  merely  to  check  the  action  of  certain  operators  on 


smooth  functions. 


3 


2.  Convergence  of  Finite  Dimensional  Distributions 


Let  9)  denote  a probability  space,  { 9^}  a nondecreasing 

sequence  of  sub  o-algebras  of  9,9  the  space  of  progressively  measurable 

real  valued  processes  f on  [0,°°),  adapted  to  { 9^)  and  such  that 

sup  E | f ( t ) | < ®.  Let  f and  f be  in  9.  Define  the  limit 

"p-lim"  by  p-lim  f = f iff  sup  sup  E | f ( t ) | < 00  and 

n n t 

Elf  (t)  - f(t)|  •>  0 for  each  t.  For  each  s > 0,  define  the 
n 

operator  9(s):9-*y'  by  9(s)£  = function  in  9 whose  value 

at  t is  the  random  variable  E _ f(t+s).  There  is  a version  which 

y t 

is  progressively  measurable  ([1],  Appendix)  and  we  always  assume  that 

this  is  the  one  which  is  used.  The  9(s)  , s > 0,  are  a semigroup  of 

linear  operators  on  9.  Let  ¥ denote  the  subspace  of  9*  of 

p-right  continuous  functions.  If  the  limit  p-lim  [-  ( s)  f — f ) ] 

s-*0  s 

exists  and  is  in  9Q,  we  call  it  Af  and  say  that  f e S>(A)  . 

The  operators  i^(s)  and  A are  analogous  to  the  semigroup  and  weak 
infinitesimal  operator  of  a Markov  process.  Among  the  properties 
to  be  used  later  is  ([1],  equation  (1.9)) 


9(a)  f - f = y(i)  Af  dT,  f e 0(A), 
0 


(2.1a) 


or,  equivalently 


From  [1],  with  slightly  altered  terminology.  Sometimes  we  write  f^ 
and  sometimes  f(t)  for  the  value  of  a process  f at  time  t. 


r 


i 

i 

i 

i 

i 

i 

t 

i 


(s 


J f(t+s)  - f(t)  » Ea  Af (t+ 1 ) d l j for  each  t ^ 0. 
rt  J o *t 


(2.1b) 


If,  for 

some 

process  Z 1 ( • ) , = 

°<v 

s < t) , we  may 

* 

write 

.—m  t!  L 

SF  , T 
t * t 

and 

Al  for  J?t,  >(  t) 

and 

A,  resp.  Let 

C and  C1 
o o 

denote  the  spaces  of  real  valued  functions  on  Rm  which  vanish  at  <*' 
and  which  are  continuous  and  which  have  continuous  partial  derivatives 
up  to  order  i (and  which  also  vanish  at  infinity) , resp.  Let  C and 
C1  denote  the  sets  of  these  functions  which  have  compact  support. 

The  following  Theorem  (a  specialization  of  [1],  Theorem  3.11)  is 
our  main  tool  for  dealing  with  (1.1).  Henceforth,  unless  otherwise 
mentioned,  t » 0 replaces  n » 00  in  p-lim. 


t 

11 

[ 

I 
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Theorem  1 . Let  Z ^ ( * ) = x L ( • ) , y t ( • ) , c > 0,  denote  a sequence 
of  Rm+m  valued  right  continuous  processes,  x(*)  a (Rm-valued) 
Markov  process  with  semigroup  Tfc  mapping  into  Cq  and  which  is 

strongly  continuous  (sup  norm)  on  CQ . For  some  X > 0 and  dense  set 
D in  c let  Range  ( A —A | D ) be  dense  in  (sup  norm, 

A = infinitesimal  operator  of  x(*)).  Suppose  that,  for  each  f t D, 
there  is  a sequence  {fc}  of  progressively  measurable  functions 
adapted  to  ( and  such  that 


p-lim[fe  - f(xC(-))J  = 0 


p-limlALfl  - Af  (x L ( • ))  ] = 0. 


(2.2) 


Then,  if  xL(0)  ■*  x(0)  weakly,  the  finite  dimensional  distributic 
of  x‘ (•)  converge  to  those  of  x(-). 


(2.3) 


Equations  (2.2)  and  (2.3)  are  equivalent  to  (the  limits  are  taken 
for  each  t as  i * 0) 


sup  E | f 1 (t)  - f ( x 1 ( t ) ) | v E | f L (t)  - f(xL(t))  I - 0 (2.2') 

«• , t 

sup  E|A!f‘(t)  - A f ( x 1 ( t ) ) I v .v,  E | A 1 f 1 ( t ) - Af(xC(t)>|  -»0.  (2. 3') 

t 


3 . Tightness 

Let  y(*),  y ‘ ( • ) , \ (•)  denote  the  functions  in  the  model  (1.1) 

or  (1.2).  Let  ^ and  denote  the  (oompleted)  o (y  ; s v t) 

and  o(y‘,  s - t).  Write  E and  El  for  E and  E -£,  resp. 

J t t 

Again,  we  describe  results  from  (1].  Let  Dm(0,"’)  denote  the 

space  of  Rm  valued  functions  on  10,"’)  which  are  right  continuous 

on  10,"')  and  have  left  hand  limits  on  (0,"’).  Note  that  x*  ( •)  e 

Dm  | 0 t kV  ) w.p.l.  Suppose  that  the  finite  dimensional  distributions 

of  x‘(*)  converge  to  those  of' a process  x(-),  where  x ( - ) has 
m 

paths  in  n 10,"')  w.p.l.  Then,  as  noted  in  111,  bottom  of  page  628, 

(x'(.)l  is  tight  in  Dm[0,"')  if  lf(xL(*))}  is  tight  in  D[0,"') 

' + ' 3 

for  each  f i C.  (C  is  used  there,  but  it  can  be  replaced  by  C .) 

It  follows  from  [11,  Theorem  4.20,  that  |f(x  (•)))  is  tight  in  D[0,"') 
it  x ' » x q weakly  and  if,  for  each  real  T ' 0,  there  is  a random 

variable  > ( ) such  that 

Et\('S)  - F‘t  min|  If  ^t+u1  " 13.1) 

for  all  OvtvT,  0 mi  \ ^ 1,  and 

+or  by  any  set  of  functions  dense  in  C in  the  sup  norm. 
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lim  lim  EY  (6)=0.  (3.2) 

5-0  €-0 


In  (1],  p.  629,  Kurtz  suggests  a method  of  getting  the  Y£(6). 

This  method  is  developed  in  Lemma  1 and  is  used  in  the  sequel.  The 

fl  below  will  be  obtained  in  the  same  manner  as  we  will  obtain  the 

of  Theorem  1.  We  have  (||f||  = sup|f(x)|) 

x 


+ 


^X+u> 


f2 (X^) 


(3.3) 


A 3 ^ 

Lemma  1 . Let  f t;  c ,and  let  there  be  a sequence  {f  } in 

where  (fL)x  £ £?(AL)  , i = 1,2,  and  such  that,  for  each  real  T > 0 

there  is  a random  variable  M such  that 


sup|fC(t)  - t (x'l  | - 0 w.p.l,  as.  e - 0 (3.4) 

t<T 


sup  lA'tf^t))1!  < M,  w.p.l,  i = 1,2. 
e>0,t<T 


(3.5) 


Then  (f(xL(*))}  is  tight  in  D[0,“). 


Proof . By  (2.1) 

(f C (t+u)  ) 1 -(fE(t))X  = jUE^A£(fe(t+T))1dX, 

from  which  (3.3)  to  (3.5)  yield  a Ye(<5)  of  the  form  (for  the  interval 
[ 0 , T ) ) max { 1 , ) where 
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& * M.  [ 6 + l sup|  (f L (t)  ) 1 - fi(xh  | , 

1 i*l  t <T 

where  is  a random  variable.  Q.E.D. 

4.  Assumptions  for  Model  (1.1);  Bounded  y(*)- 

( Al ) f ( • , • ) and  G ( • , • ) are  continuous,  and  the  first  partial 
and  first  and  second  x-partial  derivatives  of  G and  F, 
resp. , are  continuous. 

(A 2)  There  is  a constant  M sue))  that 

| F (x , y ) | + |G (x,y)  j < M ( 1+ | x | ) . 


(Al)  and  (A2)  assure  the  global  existence  of  solutions  to  (1.1). 

( A J ) y ( • ) is  a right  continuous,  bounded  stationary  process  and 
EF(x,yg)  = 0,  each  x. 

There  is  a measurable  function  P ( • ) such  that 

sup  |P(B, |B  ) - P(B-) | v p(i) , 

®i't 

where  e o(yu,u  < t),  B2  e e(yu,  u 2 t + x)  ' 

[ p(t)dt  < “.  (4.1) 

' 0 

(4.2)  is  weaker  than  the  conditions  on  the  mixing  rates  in  [4]. 

Define  the  operator  A on  C2  by  (the  subscript  x denotes  gradient) 
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Af  (X) 


KG’  ( x , y ^ ) f x ( x ) + 


vV 


dl  EF'  (x,yg)  (F'  fx(x)  )x' 


(4.2) 


I bi  (x)  f (x) 

i ' X i 


i *. 


i , i 


a . ( x 1 

l.l 


*-!!*) 

.'x  . ?x  • 
i j 


By  i A 1 ) and  ( A J ) , the  integrals  on  the  right  exist.  in  tact,  the 

fT 


improper  Lebesgue  integral  is  absolutely  convergent  ti.e.  E(  ) dt 

* ‘ , , ' 0 
converges)  unitormly  in  x . Furthermore,  if  t i c then  At  (x)  is 

continuously  d i f ferentiable  in  x,  and  the  gradient  of  At  (x)  is  the 

function  obtained  by  simply  replacing  the  argument  of  K in  1 4 . _’ ) 

by  its  x-gradient.  If  this  is  done,  then  the  improper  Lebesgue 

integral  still  is  absolutely  convergent,  uniformly  in  x. 


(A4)  A i s t he restriction  to  C*"  of  the  strong  .infinitesimal 

ope rat or  ot  a^  strong  Markov  conso r vat  ive  ( no  t i n i t e 
escape  time)  diffusion  process  x('),  with  semigroup  T . 
T(  maps  into  and  is  strongly  cont inuous  on  C ^ . 

(A5)  t ) is  right  continuous.  l.e.,  ^ * 0 ^ + ^ » 

each  t ^ 0. 

(Abl  The  set  (\-A)C2  (g:  g (\-A)f,  t t C3)  is  dense  in 

C ^ for  some  \ v 0 . 


This  follows  readily  from  the  strong  mixing.  It  is  also  a con- 
sequence of  Billingsley  [7],  p.  170.  by  using  EF(x,y,)  0,  t4.1) 

and  the  fact  that  the  functions  have  bounded  supports(let  r - 1, 
s * "•  in  (7),  equation  (20.2.1),  with  a proper  identification 
ot  s , n there).  Wo  will  use  this  and  similar  facts  frequently  in 
the  sequel. 


Remark  on  (A5) . Let  fee  and  let  f denote  either  Ff,Gf 
of  any  of  the  g ^ or  f ^ introduced  below.  Condition  (A5)  is 
introduced  only  because  we  want  E^+gf (x^+g , t+s+u)  to  converge  to 
EtT(xt't+u)  in  probability  as  s I 0.  Many  of  the  calculations  in 
Theorems  2 and  4 involve  this  type  of  right  continuity  together  with 
uniform  integrability. 

Remark  on  (A6) . Some  condition  such  as  (A6)  is  required  for 

use  of  Theorem  1.  Let  Ac  denote  the  strong  infinitesimal  operator 

of  Tt  acting  on  CQ.  Then  (A6)  is  equivalent  to  Ac  being  the 

closure  of  the  operator  A (of  (4.2)  acting  on  C , or  on  C , 

since  c is  dense  in  C:  in  the  norm  | | f | 1 2 = sup(|f(x)|  + 

x 


fx(x) 


+ |f  (x)  |))  . This  condition  does  not  seem  to  be  particularly 

XX 


restrictive.  We  only  remark  that  it  holds  in  the  two  extreme  cases: 
(1)  the  b^(’)  and  a^(-)  in  (4*2)  are  bounded,  satisfy  a uniform 
Holder  condition  and  a(*)  is  uniformly  positive  definite;  (2)  where 

Tt 

respect  to  the  norm 


2 2 2 
maps  Cq  into  Cq  and  is  strongly  continuous  on  C Q with 


| f | | 2 defined  above.  (The  same  remarks  were 

made  by  Kurtz  [1],  p.  632.) 

.can 

In  case  2,  we  /actually  consider  Tfc  as  acting  on  the  Banach 
2 

space  Cq  with  norm  | |f | \^,  and  where  f (x)  and  its  first  and 


second  derivatives  go  to  zero  as 


00,  and  modify  (A4)  accordingly. 


In  this  case,  the  closure  of  the  operator 


(the  domain  of  A is 


C here)  is  just  the  strong  infinitesimal  operator  (of  Tfc)  acting  on  its 
2 

domain  in  CQ.  Suppose  that  there  is  a matrix  valued  o(>)  such  that 
o(x)o'(x)  = a(x),  that  (A4)  holds  (as  modified  above)  and  that 
bf  ( • ) , are  locally  Lipschitz . Then  (see  remarks  below  on 

bounding) , it  is  enough  to  prove  Theorem  2 under  the  additional  condition 
that  b^(*)»  °ij(‘)  are  bounded  and,  for  arbitrary  N,  arbitrarily 


1 J 


10 


smooth  out  of  the  sphere  SN  of  center  0 and  radius  N.  Assume  that 
( • ) and  °ij(*)  have  continuous  first  and  second  derivatives.  Then 
by  the  above  remark  on  bounding,  we  can  assume  that  the  coefficients  are 
bounded  on  Rm.  Then  (Gikhman  and  Skorokhod  [13],  Chapter  8.4), 
case  (2)  above  holds.  The  conditions  imposed  are  weaker  than  those 
in  [4]  when  F and  G do  not  depend  on  t. 

Remark  on  bounding  the  coefficients.  Suppose  a(x)  = o(x)a'(x), 
where  b^  ( • ) and  ctj(*)  are  locally  Lipschitz  continuous  and  x(*), 
the  diffusion  with  generator  A,  is  conservative.  Define 

N]  M 

an  N-truncation  as  follows.  Let  b^(*),  cr  ( • ) equal  b^  ( • ) , 0—(’), 

resp. , in  S^,  be  bounded  on  Rm,  have  bounded  derivatives  of  any 

, at  least 

desired  order  in  the  complement  of  S_„  and  be/as  smooth  in  S„„,  - 

2N  2N  N 

as  b^(-),  are*  Then  the  Ito  process  xN(')  with  coefficients 

b^  ( ‘ ) , o^_.  ( • ) is  called  an  N-truncation  of  x ( ■ ) if  the  k?1  ( • ) and  oN  ( • ) 

can  be  obtained  by  a modification  of  F and  G in  Rm-Sri . N-truncations 

always  exist  since  we  can  multiply  F and  G by  suitable  smooth  real  valued 

functions  m^(x)  and  (x) , resp.,  which  equal  unity  in  S^.  We  only  need 

prove  Theorem  2 and  verify  (A4)  and  (A6)  for  sane  N-truncation  for  each  N. 

The  proof  of  the  assertion  will  be  emitted.  It  is  essentially  a note 
N 

that  the  parts  of  x (•)  and  x(")  up  to  the  first  escape  frem  are 

equal,  and  that  the  probability  of  escape  from  SN  on  an  interval  [0,T] 
goes  to  zero  as  N -+•  °°,  for  each  fixed  x^.  Here,  we  suppose  that  x(-)  are 
defined  with  respect  to  the  same  Wiener  process. 

5 . Proof  of  Weak  Convergence;  Bounded  y(»)  and  (1.1) 

The  main  job  in  using  Theorem  1 is  to  get  when  f is  given. 

To  do  this,  we  use  an  idea  exploited  for  a similar  purpose  in 
Section  3 of  [5]  and  in  Section  4 of  [3],  We  look  for  functions  of  the 
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form*  f (x,t)  = 
Define  operators 

function  of  x 
for  each  x.  At 


f (x)  + rf^U.t)  + t-  ^f  2 (x , t ) 

~ t " t 

A and  A as  follows, 
x y 

and  such  that  g(x,t)  is  a 
x = x£,  y = y^ » let 


, and  define 
Let  g(x,t) 

function  ol 


f L (x£, t)  = f 1 (t) . 
be  smooth  as  a 

V s - 


A^g(x,t)  = g^(x,t)  [ -- ( 1 + G ( x , y ) ] ; (5.1) 


i.e.,  A^  is  A , but  acting  on  g(x  ,t)  considered  only  a function 
ot  its  first  argument.  Let  A^g(x,t)  be  A g(x,t) , but  where  g 
is  considered  to  be  a function  of  its  second  argument  only.  Assume  for 
the  moment  that  A = A^  + a\  Then,  in  order  to  use  Theorem  1,  we  apply 
A to  f 1 ( x , t ) , and  insist  that 

[fx(x)  + efl,x(x,t>  + e2f2,x(x,t)1  ’ + G(x,t)  ] + 

+ Ay(ef^(x,t)  + t;2fl2(x,t))  (5.1*) 

= Af (x,t)  + 0(0, 

where  equation  ( 5 . i ’ ) must  determine  both  the  operator  A and  equations 
yielding  the  f^(x,t).  In  Theorem  2,  we  merely  write  down  formulas  for 
the  t.(x,t)  and  verify  the  conditions  of  Theorem  1. 

Theorem  2.  Under  (Al)  to  (A6) , x'(-)  converges  weakly  in 
Dm  l 0 , 00 ) to  the  di  f fusion  x ( • ) , with  initial  condition  x . 


Proof.  First  (Parts  1 to  3)  we  prove  convergence  of  finite 
dimensional  distributions,  using  Theorem  1.  In  Part  4 tightness  is 
proved,  via  Lemma  1.  Henceforth,  f is  a fixed  element  of  C 5 . 

Since  C3  is  dense  in  C2  in  the  norm  ||f|l2'  is  enough 
- 3 

to  work  with  C . 

*For  each  x and  t,  fj(x,t)  will  be  a function  of  y^,  s * t. 
The  discussion  in  this  paragraph  is  purely  formal. 


i: 


Part  1.  f^(x,t)  is  defined  to  be  a solution 

(5.1'J)  to  AyfJ_(x,t)  = -g:  (X,y^)  = -F' (x,y^)  fx(x)  . 
define  f ^ (t)  = f^(x^,t),  where 


f J(x,t) 


gm 

E^g.  (x,y  {—  + s)  )ds 

i0  1 zz 


(suggested  by 
More  precisely, 


(5.2) 


,E^l(X'yt+s)  ds 


(both  forms  will  be  used) . The  improper  Lebesgue  integral  exists  and 
is  bounded  and  absolutely  convergent,  uniformly  in  u>,x  and  in  t in 
bounded  sets,  by  the  strong  mixing  (A3),  and  the  facts  that 

EF(x,yg)  = 0 and  that  g^  has  compact  x-support.  Furthermore,  there 

e £ £ . 

are  versions  of  f^Xjt)  and  f1(xfc,t)  which  are  progressively 

measurable . 

We  next  show  that  f^(t)  t £2(AL)  . We  have 

ALfJ(t)  = P-lim  [E^(x^+6,t+6)  - f*(x*,t)]/6  (5.3) 

" p-lim  [Ef{f^(xf  - » t+6)  - fe(x^,  t+<S)  >1/6 
,5-*0  r 1 t't’0  1 '~ 

+ p-lim  [E^f L (x!;  ,t+6)  - f.L  (x^,t)  ]/6 
6-*-0  1 C 1 

if  the  limits  exist  and  are  in  & . it  is  easy  to  verify  that  the 
second  limit  exists,  is  in  &Q  and  equals  -gL (x^ ,y£) /e  . Now, 
fj(x,t)  is  differentiable  in  x.  Indeed, 


fl,x(x't)  = J2  J0Etg.l,x(x'yt+s)ds' 


since 


£ £ 

|E  g.  (y  )|ds  converges  uniformly  in  x,  and  in 

t 1 t X ltS 


t in  bounded  set*, as  T -*  “.  This  fact  together  with  the  representation 
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Et(f  1 (xt+6* t+(S)  ~ t+tf)  J/^  = 


e , £ 


(5.4) 


k |*  Etlil,»lllt«.'t*11 1 ' iFUttutY-—  * G(x;+u,y;tu)]du 

and  the  facts  that  f!”  is  right  continuous  in  the  mean  at  t,  and  that 

1 / X 

the  integrand  is  zero  out  of  a bounded  x,v  range,  can  be  used 
to  show  that  the  first  limit  on  the  right  side  of  (5.3)  exists,  is  in 
_y'0  and  equals  (x  = x^ , y = ylfc) 


If.1  ( x , t ) ] • lP(x,y)/e  + G (x,y)  ] 

■L  » X 


=j  A^  f ^ ( x , t ) . (5.5) 


Part  2.  With  x = x^,  we  will  define  f^(t)  = f2(x,t),  where 
f^(x,t)  is  the  formal  solution  to 


Ayf^U.t)  - -g2  (x , t) 

" -[F'(x,yJ)fJ>x(x,t)  + GI(x,y^)fx(x)  - Af(x)l  , 

where  A is  defined  in  (4.2).  More  precisely,  define  f2  by 


(5.6) 


f 2 (x , t)  = 


r 


Etg2 (x ' t+s) ds. 


(5.7) 


There  are  versions  of  fj(x,t)  and  fjtx^t)  which  are  progressively 
measurable.  We  now  ignore  the  G terms , for  the  difficulty  lies  with 
the  Ff,  term.  The  improper  Lebesgue  integral 

| ELg2 (x, t+s) | ds  of  (5.7)  converges  absolutely,  uniformly  in  w.x 
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and  in  t in  any  bounded  set,  by  the  strong  mixing  property  and  the 
definition  of  A.  (Indeed,  this  is  the  reason  for  the  choice  of  A.) 
To  see  this,  note  that  (5.7)  (without  the  G-terms)  equals 


0 ds{EtP'  (x'yt+s)fl,x(x't+s)  - EF”  (x'yt+s)fi,x(x't+s)  (5-8) 


where  the  average  value  of  the  coefficient  of  ds  is  zero  (by 

stationarity  and  the  definition  of  A),  and  use  the  strong  mixing 

2 

condition.  In  fact,  using  the  change  of  variables  s/ 1 “ *■  s',  it 
can  be  seen  from  (5.8)  and  the  strong  mixing  and  compact  x-support 

of  f^  that  | f ^ ( t ) | is  bounded  w.p.l,  uniformly  in  x and  u> 
in  bounded  t intervals.  Furthermore,  g,(x,t)  is  continuously 
differentiable  in  x.  In  fact,  the  convergence  assertion  in  the 
sentence  above  (5.4)  also  holds  for  g^  replacing  g^  and 


and 


I2,x(x,t)  t 2 


Etg2  x.(x't  + s)ds- 


Expression  (5.3)  holds  with  f2  replacing  f^  if  the  limits 

A 

exist  and  are  in  , Again,  we  readily  verify  that  the  second 

limit  in  (5.3)  exists,  is  in  ^ and  equals  - g2(x,t)/t"  (x  - xfc) 
An  argument  similar  to  that  leading  to  (5.5)  yields  that  the  first 
limit  in  (5.3)  exists,  is  in  and  equals  (5.5)  with 

replacing  f^. 

Part  3.  Now,  we  apply  Theorem  1.  Since 


sup  E(|fJ(t)  | + | f i (t)  | ) < », 

t,t^o  1 * 
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we  have 

p-lirnff1-  f(xL(’))J  = 0. 

A f ^ t t 

Now,  calculate  A f . By  Parts  1 and  2,  with  x = xfc  , y = yfc, 

Afl(x,t)  * A Lf ( x ) + £Alf^(x,t)  + t^A^f ^ (x , t)  (5.9) 

= [ F ( x , y ) / 1 + G(x,y) ] ' £x (x) 

+ el-F'(x,y)f(x)/t2  + (F(x,y)/t  + G(x,y))'fl  (x,t)] 

X 1 i X 

+ t2(-  {F ' (x , y ) f j x<x,t)  + G ' (x , y) f x (x)  - Af (x) } 

+ (F(x,y)/e  + G<x,y))'f*  v(x,t)))  = 

^ » X 

= Af (x)  + t ( G ' ( x , y ) f ^ x(x,t)  + F ' (x,y) f ^ x(x,t)] 

+ t2G' (x , y ) ( x , t ) . 

We  now  readily  verify  that  p-limlA^f^  - Af(xl(*))l  = 0.  Since 

x0  = Xq,  Theorem  1 implies  that  the  finite  dimensional  distributions 

converge . 

A 

Fart  4.  Tightness.  For  tightness,  we  use  Lemma  1.  Each  f c C 

* 3 

can  be  approximated  uniformly  arbitrarily  closely  by  an  f t C . Thus, 

by  the  Lemma  and  discussion  preceeding  it,  we  only  need  prove  that 

. /\  3 \ 3 

{ f ( x (•))}  is  tight  for  each  f e C . Let  f t C and  construct 

f^,f‘  exactly  as  the  f^f^  were  constructed  above,  and  define 

f 1 ( t ) as  above.  Then  fL  and  (fL)2  are  in  £j*(Ak)  and  there  is 

a constant  M such  that  w.p.l 
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sup  i = 1,2, 

e>0,w,  t 

sup|ft(t)  - f (xL  (t)  ) | + 0,  as  e -*■  0. 
t 


Thus,  Lemma  1 implies  that  {f(xL(*))}  is  tight  in 

~3  c > m 

each  ft  c ; hence  {x  (•)/  is  tight  in  D [0,°°). 

together  with  the  convergence  of  finite  dimensional 

implies  weak  convergence.  Q.E.D. 


D [ 0 ,°°)  , for 
The  tightness, 
distributions 


An  approximation  to  an  integral.  In  problems  where  changes  of 
measure  via  Girsanov-like  transformations  are  involved,  such  as  occur 


in  some  identification  problems  (Balakrishnan  [ 8 ] , [ 9 ] ) , we  need  to 

rt 


get  limits  of  integrals  such  as 
ii 


q ’ (Xg)(Yg/e)  ds , as  t -*  0. 


Let 


denote  the  real  valued  functions  on 


-m+m’+l 

R with  compact 


support,  whose  i 

< - ;t 


, th 


(y  /t)ds  and 


partial  derivatives  are  continuous.  Let 

t. 


e , e c t. 

u = (x  , Y , z ) . 


Then  ut  is 


„m+m ’ +1  , , , 

R valued  and 


uL  = G(xt,yt)  + i F(xL,yt), 


where 


G(x,y)  = (G (x , y) ,0,0) 

F(x,y)  = (F(x,y) ,y,q’  (x)y)  . 


The  remarks  below  (A6)  all  pertain  here  also. 
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Theorem  j . Let  q(-)  satisfy  the  conditions  on  F(*)  ui 

( A 1 ) — ( A2 ) . Assume  (Al)  to  (A6)  where  (A4)  and  (A6)  hold  for  the 

— '2 

process  u(*)  and  operator  A defined  on  C by 


Af(u)  = EG’ U,y0) fu(u)  + j dlEF’ (x,y0) (F* (x,yt) fu(u) )u.  (5.10) 


Then  u'' ( • ) converges  weakly  in  Dm+m  ^lO,**)  to  u(*)  = (x(’)r 

\ 2 

Y ( • ) , z ( • ) ) » uq  = (x0»0f 0) , a diffusion  with  generator  A on  C . 

r 

The  process  Yfc  is  a Brownian  motion  with  covariance  K - 2 Ey  y\U. 

/Ito  ° 

The  limit  z ( • ) has  the  representation  (the  expectation  is  over  y ( • ) ) 


in  terms  of  the  limits  x(')#Y(‘) 


dz  = q’(x)dY  + | EF  ’ (x  , yQ )(q  1 (x)  y t ) ^d  l ] dt . 


(5.11) 


(The  last  term  in  (5.11)  is  the  so-called  correction  term  of  the 
limiting  integral  approximation.) 

The  proof  of  weak  convergence  of  ul(*)  to  u(-)  is  simply 
an  application  of  Theorem  2.  Once  the  weak  convergence  is  known, 
then  the  representation  (5.11)  is  not  hard  to  get,  and  we  omit  the 
details. 


6.  Unbounded  y ( • ) and  (1.1). 


Our  approach  here  will  be  only  a little  different  from  that  in 
Section  5.  In  order  to  avoid  conditions  which  look  overly  complicated. 


we  specialize  F(x,y)  and  G(x,y)  to  F(x)y  and  G(x),  resp. 


1 8 


Assumptions.  In  this  section,  the  convergence  of  finite  dimensional 
distribution  is  proved,  and  tightness  is  treated  in  the  next  section, 
owing  to  the  unboundedness  of  y{*)»  it  is  convenient  to  artificially 
bound  the  F and  G . We  do  this  by  dealing  with  a sequence  ot 
approximations  to  the  original  processes.  The  operator  A is  still 
defined  by  (4.2). 

(Hi)  F ( • ) and  its  first  and  second  order  partial  derivatives 
are  continuous. 

( H 2 ) G ( • ) and  its  first  order  partial  derivatives  are 
cont inuous . 

IB3)  1 ) (again  completed)  is  right  continuous,  ami  so 

is  the  stationary  process  y(*),  w.p.l.  (See  the  remark 
concerning  ( AS ) . ) 

pef ine 


vt  " t?  j0<Wis'  vt  - |0Etyt+sds- 


(H4) 

For 

some  P s 0 , sup  E ( j 

' olEtyt+Jds) 

Thus , 

and  v are  well-defined. 

(B5) 

Oyt 

2+p 

| < »,  some  p N 

0. 

(B6) 

sup 

t 

E yt+sVt>sl) 

Note 

that 

"'Wt*.  ■ J 

f 0Et+s^  t + s+u>' 

2+P 


V >»•  , 


2+p 


v <«•,  smie  *'  s 0. 


" f0Eytyi*udu  *nd 


does  not  depend  on  t or  s,  and  is  well-defined  by  (H4),  (H5) . 
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Owing  to  the  special  form  of  F and  G,  there  are  locally 

Lipschitz  b(*)  and  o(‘)  such  that  (see  (4.2))  o(*)o'(*)  = a(*). 

N C N 

See  also  the  remarks  on  bounding  below  (A6) . x and  x ' (•) 
denote  the  N-truncations  of  x(-)  and  xL(*)- 

(B7)  For  a sequence  N 00 , there  are  N-truncations  which 
satisfy  (A4)  and  (A6). 


Remark  on  assumptions  (B3)-(B6).  Let  w(’)  denote  a vector 
Brownian  motion,  Q a matrix  with  eigenvalues  in  the  open  left  half 
plane  and  let  D and  G be  matrices.  Define  processes  Y(*)  and 
y ( • ) by 


dY  = QY  dt  + Ddw, 
y ( • ) = GY  ( • ) . 


(6.1) 


Then  (B3)  to  (B6)  hold.  In  this  case,  we  can  let  and  & 

2 

--J  Y , s < t,  resp. , in  all  the  foregoing, 
s - 


measure  Y , s < t/e  , and 
s - 


Then  v is  proportional  to  Yfc  and 
c^  are  positive  constants. 


c s 

lEtyt+s'  - I Yt ' cle  2 ' where  the 


Theorem  4 deals  with  the  convergence  of  finite  dimensional 
distributions. 

Theorem.  4 . Under  (Bl)  to  (B7)  and  the  first  sentence  of  (A4), 

£ 

the  finite  dimensional  distributions  of  x (•)  converge  to  those 
of  x ( • ) with  initial  condition  xQ,  as  e 0. 
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6.  N 

Proof.  If  the  finite  dimensional  distributions  of  x (•) 
converge  to  those  of  x^  ( • ) (initial  condition  x^)  as  e ^ 0 for 

a sequence  N -*•  00 , then  the  conservative  and  the  strong  Markov 
properties  of  (A4) , (B7)  yield  the  theorem.  So  we  only  need  prove 


convergence  for  a fixed  N.  Consequently,  we  may  assume  that  F 
and  G are  bounded  and  drop  the  affixes  N. 

The  details  are  very  similar  to  those  of  Theorem  2, 

€ £ 

and  we  only  make  a few  remarks.  As  before,  define  ^ and  f ^ by 
f^(t)  = fj-(x£,t),  f^(t)  = f^Cx^.t), 

t 1 

where  f^(x,t)  is  defined  as  in  Theorem  2.  These  functions  are 

£ 

no  longer  bounded,  but  still  sup  E | f . (t) | < ».  It  is  rather 

t,e>0  1 

straightforward  to  verify  (in  fact,  easier  than  in  Theorem  2 owing 
to  the  special  form  of  F(x,y)  and  G(x,y)  here)  via  (Bl)  to  (B6) , 
that  ff  £ 9(^1  and*  take  the  same  values  as  in  Theorem  2. 

Furthermore,  the  expectations  of  the  absolute  values  of  the  coefficients 
of  £ and  £2  on  the  far  right  side  of  (5.9)  are  bounded,  uniformly 
in  t.  Also 

sup  E{  | A^f  (t)  | + | f ^ ( t ) | } < <*>. 

’ t,e>0 


Thus , 

* ' — — — 

For  example,  to  show  that  the  expression  for  *efF(t)  is  in  ^ we 

note  that  the  compact  support  of  f and  (B4)  to  (B6)  imply  uniform 

integrability  of  the  expression.  This,  together  with  (B3)  yields 

p-right  continuity. 
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p-lim[fL  - f (x1!1)) 1 = 0 and  p-lim[Atfc  - Af(xL(*))l  = 0, 


from  which  the  theorem  follows,  by  Theorem  1.  Q.E.D. 


7.  Tightness;  Unbounded  v(') 


Owing  to  the  unboundedness,  it  is  more  difficult  to  prove 
tightness  via  the  method  of  Lemma  1.  To  avoid  (what  are  at  the 
moment)  awkward  conditions,  we  suppose  that  y(>)  satisfies  (6.1). 
Then  the  f*-  can  be  explicitly  evaluated  and  the  proof  is  easy. 

and  that  y ( • ) 

Theorem  5.  Assume  (Bl) , (B2) , (B7),  the  first  sentence  of  (A4)  / 
satisfies  (6.1) . Then  { x C ( • ) } is  tight  and  converges  weakly  to 
x ( • ) as  e •»  0 . 


Remark . The  tightness  argument  only  uses  the  compact  support 
of  f , (Bl) - (B2)  and  (6.1). 


Proof . The  method  and  notation  of  Lemma  1 and  Theorem  2,  part  4, 
.11  be  used  here.  We  need  to  show  that,  for  each  T and  each 


f t C3 


(i)  lim  sup|f(x.)  - f (t) | = 0 w.p.l 
£-*0  tvT 


(ii)  (fL)2  t ^(A^) 


(iii)  lim  sup  | AL  (f  **  (t)  ) | < °°,  j = 1,2,  w.p.l. 

£>0  t<T 


In  our  case,  there  is  a matrix  Cn  such  that 
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fj(x,t)  = (F(x)C0Y(t/£2) ] ’fx(x) 

g2(x,t)  = {F(x)y(t/e2) } ' { [F  (x) CQY (t/t 2) ) ' fx (x) )x 

+ G ' (x)  f (x)  - Af (x)  . 

X 

f^(x,t)  is  a quadratic  form  in  the  components  of  Y(t/t  ) where 
the  coefficients  are  bounded  differentiable  functions  of  x with 

C 2 A £ 

compact  support.  Also,  as  is  readily  verifiable,  (f  (t)  ) £ £#(A  ), 

and  AL(f L (t) )3 , j = 1,2,  have  terms  in  powers  of  the  components  of 
2 

Y (t/£  ) up  to  2j  + l.  Thus,  to  verify  (i)  and  (iii)  it  is  enough 
to  verify  that  for  each  T > 0 

lim  sup  e|Y(t/t2)|3  = 0 w.p.l.  (7.1) 

e-*0  t<T 


Equation  (7.1)  holds  since,  for  each  a > 0,  the  Gaussianess, 
stat ionar ity , and  special  form  (6.1)  imply  that  there  are  finite 
w.p.l  w-functions  and  such  that  | Y ( t ) | < C^t  + C2  for 

all  t,  w.p.l.  Q.E.D. 

8 . An  Approximate  Jump  Case;  (1.2)  With  No  yL(*)  Term 

Since  the  classical  papers  of  Wong  and  Zakai  [6) , the  problem 
of  using  Ito  or  other  types  of  equations  to  approximately  model 
processes  which  are  the  solutions  to  ordinary  differential  equations  has 
received  much  attention;  e.g.,  ( 3 1 — [ 5 ] and  Sections  5 to  7 above. 
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Alternative  approaches  have  been  taken  by  McShane  (10]  and  Sussmann 
111]  who  sought  either  a theory  of  integration  or  a differential 
equation  and  a topology  on  the  input  functions  so  that  the  output 
is  a continuous  function  of  the  input.  The  differential  equations 
which  were  of  the  so-called  Stratonovich  form,  which,  in  fact,  are 
precisely  Ito  equations  with  appropriate  dynamic  terms. 

Little  has  been  done  when  the  input  is  an  approximation  to  an 
impulse  (its  integral  is  an  approximation  to  a pure  jump  process  ) . 
Marcus  I 12]  has  done  some  work  along  McShane' s "belated  integral" 
point  of  view.  This  work  [12]  has  some  interesting  aspects,  but 
also  a number  of  shortcomings.  The  dynamics  are  rather  special, 
being  (in  part)  analytic  functions.  This  is  a disadvantage  in  any 
approximation  theory,  where  robustness  is  a key  word.  Some  heavy 
lie  algebra  machinery  was  used,  and  the  form  of  the  results  tended 
to  obscure  the  basic  simplicity  of  the  problem.  Also,  a very 
particular  impulse  approximation  was  used  (piecewise  constant) . In 
this  section,  we  take  a simple  minded  but  inherently  natural  and  robust 
approach,  using  pathwise  approximations  and  limits.  The  limits  being 
either  ordinary  or  stochastic  differential  equations  with  impulsive 
or  jump  inputs. 

Let  N.  (dsxdy) , i = l,...,k,  denote  a sequence  of  scalar 

(re- 
valued random  measures  and  define  N.(t)  = a N.(ds*da),  where 

1 0 

N ^ ( • ) is  taken  to  be  right  continuous.  The  range  of  the  jumps  of 
N (•)  is  a bounded  set  . Each  ( • ) is  assumed  to  have  a 

finite  number  of  jumps  on  each  bounded  interval  w.p.l,  and  the 
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probability  is  zero  that  different  NL  (•}  have  simultaneous 
jumps . 

In  this  section,  we  deal  with  the  equation 

e e k e e 

x.  = G (x  ) + l h. (x  )J  (t)  , (8.1) 

r i=l  1 r 1 

£ • 
where  the  input  (t)  is  an  "approximation"  to  the  impulse  (t) . 

Figure  1 illustrates  some  ways  in  which  an  actual  integrated  input 

j? (s) ds  might  approximate  an  ideal  integrated  impulsive  input  N^t). 
0 1 

In  the  figure  a jump  of  Y occurs  at  t = tQ.  With  approximation  (1), 

J^(t)  = Y/e  on  [tQ,t0+e].  Define  K^(x,y)  to  be  such  that 

x + K^(x,Y)  is  the  solution  to  x = H^(x)Y  at  t = 1,  with  x^  = x. 

It  is  convenient  to  start  with  the  vector  of  ideal  integrated 

inputs  N ( • ) = {N^(*),  i < k}  and  to  get  the  actual  inputs 

from  this,  as  indicated  in  Figure  1.  When  the  parameter  is  e,  we  work 

-times 

only  with  paths  for  which  the  inter  jump^  of  the  vector  { ( • ) , 

i = l,...,k}  are  > e.  Obviously,  this  involves  neglecting  a set  of 

paths  whose  probability  goes  to  zero  as  e -*•  0,  and  it  has  no  effect 

on  the  limiting  process.  Thus,  for  our  "limit"  results,  only  the  case 

of  one  input  need  be  treated,  and  unless  noted  otherwise  set  k = 1 

and  drop  the  indices  i on  J,N,H  and  K.  Let  N(*)  jump  Yj  at 

t = a . , its  jth  jump  time.  Define  p^  = jf/Y.  on  [o. ,o.  + e],  and 

equal  to  zero  out  of  u [o  a.+  e]  and  define  pt  = f^peds.  Thus 

Y 3 3 t I <rs 


1. 


Now  (8.1)  can  be  rewritten  in  the  form 


*t  = G(xt)  + H(xt)YjPt  on  l°j'0j+Ll  ' (8'2) 

= c.(x£)  otherwise. 

The  value  of  pL(-)  can  depend  on  the  jump  time  and  size  and  on  the 
state  prior  to  the  jump  (and  on  the  index  i in  case  (8.1)). 


Assumptions . 

(Cl)  The  C.  and  IP  are  continuous. 

(C2)  There  is  one  and  only  one  solution  to  x = G(x) , for 
each  xQ  in  R1";  for  each  T < • , the  solution  is 
bounded  on  ]0,T]  uniformly  in  bounded  xQ  sets. 

(C3)  For  each  i,  there  is  one  and  only  one  solution  on  [0,1] 
to  w = H i (w ) Y for  each  Y t,  and  Wq  e R . This 

solution  is  bounded,  uniformly  on  bounded  (wQ,Y)  sets. 

(C4)  There  are  real  numbers  m s 0,  M(  v such  that 

mt  v pl (t)  v on  the  [o^,Oj+t)  intervals  and  pL(") 

is  continuous . 

( C5 ) I,et  0 v pl(t),  where  / C)  is  bounded  and  continuous 

on  10,1]  and  uL(s)ds  v u . Define  wl(-)  by 
> 0 

wL  = G (w) M 1 ( t ) + H.(w)Y.  For  each  i,  let  w‘ ( • ) exist 
and  be  bounded  on  10,1]  uniformly  on  bounded  (Y,w^) 
sets  and  in  {t,u  (•),  L *■ 


t } f or  some  >-  ^ 


0. 
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Hemarka.  (C2)  implies  that  x ( * ) is  continuously  dependent  on 
x^,  and  (C3)  implies  that  (for  each  i)  w(’)  is  continuously  dependent 
on  w^.Y.  (C2)  , (C3 ) and  (C5)  are  partially  redundant,  but  it  seemed 

easier  to  state  the  assumptions  in  this  way. 

Let  w « (w) Y and  w^  ♦ w0-  Then  (C 5)  and  (CJ)  imply  that 
|wL(t)  - w(t)|  + 0 as  c ► 0 uniformly  on  10,1)  and  on  bounded 
( Y , Wq)  sets,  and  in  U»'(*)). 


Theorem  6 . Assume  (Cl)  to  (C  5).  Let  x ( • ) be  defined  by 
ft 


Xt  “ X0  + 


G(x  )ds  + 

N 11 


l j K. 

i J R, 


(x  _ ,«) N . (du*  ds)  . 
s 1 


(8.  J) 


Let  Xy  xQ.  Then  for  each  T < *“ , 


sup  |x‘-x  | ->0  as  t 0,  w.p.l, 
ttT  1 t 

where  T *>  (0,T|  ~U  (o.  , o,  +e]  . 

1 J 3 3 

If  the  dimension  m ■ 1,  then  x c x(*)  in  the  Skorokhod 

topology  117),  Section  14). 


Remarks.  (8.3)  is  the  correct  limit  equation  - the  analog  of 
the  Wong-'Zakai  or  Stratonovich  equation  for  the  modelling  of  the 
output  of  a system  with  approximate  jump  inputs.  We  do  not  quite 
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have  weak  convergence  or  convergence  in  the  Skorokhod  metric  for 
almost  all  paths,  unless  the  components  of  the  solution  to 
w = H ^ ( w ) Y are  monotonic  in  [0,1]  for  each  Y and  w^.  To  see  why, 
let  m > 2.  Then  it  is  possible  that  , the  first  (or  any  other) 

component,  behaves  as  in  Figure  2;  i.e.,  it  does  not  approach  its 
value  at  time  1 monotonically . Let  N^*)  jump  Y at  time  o, 
and  suppose  that  there  are  no  other  jumps  on  the  interval  lo-e,o+e]. 
For  small  t,  the  first  component  of  the  solution  to  (8.1)  on 
lo,o+t]  is  essentially  a compression  of  the  [0,1]  segment  in 
Figure  2 to  the  [0,c]  interval,  together  with  a shift  of  o in  the 
time  origin.  Owing  to  the  non-monotonicity,  we  cannot  have  weak 
convergence  - or  pathwise  convergence  in  the  Skorokhod  topology,  since 
in  the  limit,  as  e -»  0,  the  compression  of  the  curve  in  Figure  2 gives 
a triple  jump  at  time  0.  Because  of  this,  the  assertions  of 
Theorems  6 and  7 are  a little  awkward. 

A main  feature  of  the  theorem  is  the  robustness  of  the  result; 
the  limit  does  not  depend  on  the  precise  form  of  the  approximations 
j'(-).  Obviously,  the  interpolation  need  not  be  over  only  an 
t-interval . 


Proof . We  need  treat  only  one  jump  and  one  H.  term,  owing 
to  the  assumptions  on  N^(*)  and  on  the  continuity  with  respect  to 
parameters  implied  by  (C2)  , ( C 3 ) and  (C5)‘.  So  return  to  (8.2)  with 
o ■;  set  equal  to  zero. 


We  change  the  time  scale.  Define  a monotone  increasing  function 

rt  . 


on  [0,t],  i (t)  * i by  di/dt 


Pt  or  x(t) 


0 


p ds . Thus , i ( t ) 
s 


1, 
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and  the  inverse  t(t  ) 


exists  by  (C4 ) . 


Def  int* 


zl  ( D 


x 1 ( t ( i ) ) • 


Then 


dzL  ( i ) 
“cTI 


G(zt(i))ut'(i)  + H(zL  ( i)  )Y, 


l < 1, 


(8.4) 


where 


U ( i) 


lpK(t( D ) J 


dt(i) 

"Ti  ' 


z 1 (0)  = xl  (0)  . 


Now  p (Odi  - t and  U (•)  satisfies  the  conditions  in  (C5)  . Let 
> 0 

xL  (0)  ■>  x(0)  as  t ■+  0.  Then  xL(e)  ■*  K(x(0),Y)  + x(0)  as  i -»  0. 
This,  together  with  the  continuity  conditions  (C2),  (C3) , (C 5)  and 
a concatenation  of  the  argument,  implies  the  theorem.  Q.E.D. 


9 . The  Jump-Diffusion  Case  (1.2);  Bounded  y(») 

We  now  return  to  the  full  model  (1.2).  Owing  to  the  possible 
non-weak  convergence  of  xl (■)  » x(*)  in  the  pure  jump  case,  due  to 

the  convergence  problem  at  the  jump  times  which  was  illustrated 

in  connection  with  Figure  2,  the  combined  jump-diffusion  case  will  be 
treated  by  a *piecin9  together ’ argument . Here  N^(-),  i * l,...,k, 
are  independent  Poisson  random  measures  with  rates  X s 0,  jump 
distributions  D ^ ( * ) (with  bounded  support  R.)  and  are  independent 
of  y(-).  With  A defined  by  (4.2),  define  the  operator  A7  on 


by 


A f(x)  - Af(x)  + l X.  If  (x+K^x.u))  - f(x)lDi(d«).  (9.1)) 
J i 


k(-)  denote  the  jump-diffusion  process  whose  infinitesimal 


operator  on 


t he  | ump 


on  C2  is  A . Except  for  the  non-convergence  problem  at 
J 

points,  xLC)  will  essentially  converge  weakly  to  x(-). 


Set  oQ  = 0 and  o = qth  jump  of  the  vector  valued  process 

N(-)  - IN  ( * ) , i v k).  Let  the  qt^1  jump  be  a jump  of  N,  (•)  and  have 
1 1 1 

value  V . Define  xl  to  be  the  solution  of  (l.l)  with  initial 

q 0 

condition  x^  and  let  x^ ( • ) (q  ^ 1)  be  the  solution  to  (1.1)  with  initial 
condition  xl (o  + t),  and  where  y'(u  +t+»)  is  used  in  lieu  of  yL(*). 

‘1  q 

Let  x (•)  denote  the  diffusion  of  Sections  4 and  5,  with  initial 

q 

condition  x,.  if  q = 0,  and  x (0)  = x . (o  -o  , ) + 

0 q q-1  q q-1 

t K (x  , (c  -o  . ) , Y ) in  general.  We  will  need  either  (Dl)  or  (D2) 

q-i  q q-i  q 

q - 1 

to  replace  (C5)  . Let  A^(Y)  denote  the  operator  on  C functions 
which  is  defined  by  (H.(w)Y)'  ^ • 

IP l ) u ' l s ) / t is  bounded  in  s v 1 and  i . The  functions  P , 

C.  and  H.  are  bounded. 

l 

(02)  u‘(s)/r  is  bounded  in  s < 1 and  t.  Each  A^(Y)(Y  t R.) 

conservat  ive 

1 s the  strong  infinitesimal  operator  of  a ( d eg  She  l^aT  el  ” 

Maikov  semigroup  mapping  (and  strongly  continuous  on)  C() 

v 2 

into  Cg.  Also  (\-A^(Y))C  is  dense  in  Cq  for  each 
Y t R^  and  i,  and  some  X * 0 (which  can  depeixi  on  Y and  i). 

Theorem  7.  Assume  (Al)  to  (A6) , (Cl),  (C4)  and  either  (Dl)  and 


(Cl)  ojr  ( D2 ) . Then  for  each  N,  ( x ^ ( • ) , q v N)  converges  to 


r 


{ ( • ) / q < N 1 weakly  in  D ( 0 , 00 ) . 


Note . The  remarks  below  (A6)  apply  here  also. 


Proof . Owing  to  the  independence  of  N(')  and  y(')  and  right 
continuity  of  y(-),  y(Oj+c+*)  has  the  properties  of  y ( - ) . Due  to 
this  independence,  the  independent  increments  property  of  N(-)  and 
the  uniqueness  and  strong  Markov  property  of  the  x(‘)  of  Sections  4 
and  5 we  can  use  a "piecing  together"  method  based  on  the  following 
assertion;  Let  a component  N (•)  jump  Y at  t = o,  with  no  other 

jumps  on  [o-e,cH-e],  and  let  xb(o)  = xt(o)  -*■  x(o)  weakly  as  e ► 0 and  define 
x (•)  for  t e (o,o+t)  by 


xs  = G(xs'ys}  + r F(xs'ys)  + \(xs)psY- 


(9.2) 


Then  (to  be  proved)  x (o+t.)  converges  weakly  to  x(o)  + K.  (x(o),Y), 

ll 

as  t ‘0.  We  will  prove  only  the  assertion,  and  the  proof  uses  a 
combination  of  the  ideas  in  Theorems  2 and  6.  For  notational  simplicity, 
let  o=0  and  drop  the  index  5.^. 

As  in  Theorem  6,  change  the  time  scale  by  defining  i(')  on 
(0,c)  and  wL(*)  by  dt(t)/dt  = pL(t),  wlO)  = xL(t(i)),  where  t(*) 
is  the  inverse  of  £(')•  Then 


dwL ( i) 


= G(wl  (i)  ,yt'(t(i)))u*'(i)  + — t—  F(wt(D  ,yt(t(0)) 


+ H(wL(i))Y,  wl(0)  = xu(0) , 


(9.3) 
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where  vi£(i)  = [pL  (t  (i  ) ) ] ~1  = dt(i)/dT.  We  need  only  show  that, 

weakly  .. 

for  fixed  Y,  w (1)  converges/to  x(0)  + K(x(0),Y)  , the 
value  at  t = 1 of  the  solution  w(')  to  w = H(w)Y,  w(0)  = x(0), 
o=0-  We  will  actually  prove  the  stronger  result  of  weak  convergence 
of  w^*)  to  w(')  in  Dm[0,l],  for  each  fixed  non-random  Y. 

First,  the  proof  under  (D2)  will  be  given.  Using  the  method  of 
Theorem  2,  let  f £ C2  and  set  U *" ( t ) = 0 and  t(i)  = t(l)  = £ for 
i _>  1 and  define  f^lw,i)  by  (t  < 1) 

.CO 

f^(w,T)  = ^ e^(t)F’  (w,yL(t(t  + s)  ) )ut'(T+s)fw(w)ds. 


e > 0 

A t used  in 

In  the  definition  of  T^  and  A',  use  Ft-Tj  just  as  FT  was  / Sections  3 to  5 . 

Set  f£(x)  = f (w‘'(i)  ) + £f  ^ (wE  (T)  , T)  . Then,  it  can  readily  be  shown 

£ £ ~ £ 

that  f (w  ( i ) ) and  f^(w  (i),t)  are  in  the  domain  of  A and  that 


A^t)  = f ’ (wt(i)  ) [G(wt(i)  ,ye(t(i)  ) )Mfc(i) 


+ F(w£(t)  ,y£(t(T)))  H...L1!  + H(wt'(i))Y] 

- F‘  (wt(T)  ,ye(t(t))f  (WC(I))  + £fi’  (wL(t)  ,T)wfc(t)  , 

W -L  f w 


which  equals  A(Y)f(w  (l))  + 0(£)  (we  dropped  the  subscript  on  A(Y)) 

This  yields  convergence  of  the  finite  dimensional  distributions  of 

wc (•)  to  those  of  w(*),  as  in  Theorem  2.  Tightness  is  also  proved 

in  the  same  way  as  done  in  Theorem  2,  completing  the  proof  under  (D2) . 

/Which  we  used  in  Theorem  2, 

Note  that  the  function  fj,  is  not  needed  here. 

Now,  we  prove  the  assertion  under  (Dl)  and  (C3). 

In  this  case  |w  (t)  | is  bounded  uniformly  in  £ > 0 and 
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t < 1 and  wL(0)  -*■  x(0)  weakly.  Thus  {wL(*))  is  tight  in 

n\  , . 

C 10,1]  and  so  is  the  function  with  values 


u (X) 

[— — F(w  ( 0 ,y  (t(t) ) ) + G(w  (x),y  (t(x)))p  (x)]dx. 


Consequently,  drawing  a convergent  subsequence  and  indexing  it  also 
by  e we  have  that  wL(*)  converges  weakly  to  a continuous 
process  w(-)»  By  using  a Skorokhod  imbedding  technique,  we  can 
assume  for  our  purposes  that  the  convergence  is  w.p.l,  uniformly 
on  bounded  intervals  and  write 


w(v) 


x (0) 


+ 


fv  - 

H (w (s) ) Yds  + lim 
0 e 


,v  c 
u ( X ) 

Jo  e 


F(we(x) ,yL(t(x) ) )dx. 


(9.4) 


We  wish  to  show  that  the  limit  in  (9.4)  iz  zero  w.p.l  for  each  v. 
If  so,  then  since  it  is  continuous  w.p.l,  it  must  be  identical 
zero  w.p.l.  Then,  by  the  uniqueness  (C3) , w(*)  = w(*)  and  the 
proof  will  be  concluded. 

The  limit  equals 


lim 

e 


V t . 

M y’  F(w(X)  , y (t(x)  ) )dX 


by  the  continuity  of  F and  boundedness  of  y(*)-  Let  a > 0. 


define  w (t)  = w(mu)  on  [ma,ma+u)  for  each  integer  m.  The 

fv  uc  CD 


F(w°t  (t)  ,yL(t(i)  ) ) dt 


difference  between  the  last  limit  and  lim 

£ 

goes  to  zero  as  a -»•  o.  Thus,  we  need  only  show  that 

/•ma+a  e _ £ /f°r  each  a. 

lim  **  - 1 F (w(ma)  ,yL  (t  ( x)  ) ) dx  is  zero  w.p.l  Now,  by  changing 

t J ma  e 


the  time  scale  back  to  the  original  one,  the  last  limit  equals 
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lim 

e 


ft(ma+a)«,#-.  ...  e 


t(ma) 


F(w(mcO_,y  (u)  ) du>  Since  t(not)  0 for  each  n as 


c -*•  0 , the  results  of  Theorem  2 imply  that  this  last  limit  is 
zero  w.p.l.  Q.E.D. 
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